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Abstract — A structure formed by combined lattices of infinitely 
long wires and split-ring resonators is studied. A dispersion equa- 
tion is derived and then used to calculate the effective permittivity 
and permeability in the frequency band where the lattice can be 
homogenized. The backward wave region in which both the effec- 
tive permittivity and permeability are negative is analyzed. Some 
open and controversial questions are discussed. It is shown that 
previous experimental results confirming the existence of back- 
ward waves in such a structure can be in deed explained in terms of 
negative material parameters. However, these parameters are not 
quasi-static and thus the known analytical formulas for the effec- 
tive material parameters of this structure, which have been widely 
used and discussed in the literature, were not correct, and it was 
the reason of some objections to the authors of that experiment. 



I. Introduction 

Meta-materials with negative permittivity and negative per- 
meability, which were first suggested in p], have attracted much 
attention recently. A meta-material which has simultaneously 
negative permittivity e and negative permeability /i within a 
certain frequency band at microwave frequencies has been in- 
troduced recently [0]. This structure consists of two lattices: a 
lattice of infinitely long parallel wires, and a lattice of relatively 
small (compared to the wavelength A in the host medium) par- 
ticles which are called split-ring resonators (SRR:s). In [Q] and 
two analytical models of SRR (similar to each other) have 
been developed for the resonant permeability at microwave fre- 
quencies. Lattices of wires at low frequencies (when the lattice 
period d is smaller than A/2) were considered as homogenuous 
dielectric media long time ago in [^] and were studied again 
recently in the low-frequency region [B], P]. At these low fre- 
quencies the negative permittivity is due to the lattice of wires 
according to the models of [g], [Q], |Q] (for waves propagat- 
ing normally to the wires with the electric field polarized along 
these wires). These results were combined in [g], [jjo[], [11|, 
[ [l2| l to form a simple model predicting simultaneously negative 
e and [i within the resonant band of a SRR. In [ |io| ] this predic- 
tion has been qualitatively confirmed by numerical simulations 
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using the MAFIA code. Dispersion curves obtained numeri- 
cally contain the pass-band within the SRR resonant band (due 
to the presence of the SRR lattice). This pass-band can also 
be predicted by the analytical model. However, the numerical 
dispersion data obtained in [ |lO| ] have not been used to extract 
the material parameters. The experimental observation of the 
negative refraction of a wave in such a structure is reported in 
[|||. The phenomenon of the negative refraction was predicted 
in |[l|] for media with e < and /1 < 0, and according to this 
theory they correspond to the backward wave region (where the 
Poynting vector of the eigenwave is opposite to the wave vec- 
tor). 

Does the experimental observation of [||] mean that the struc- 
ture suggested in [[|] can be described through e and [i which 
are both negative within the SRR resonant band? Based on ^ 
one can only assert that the backward wave region necessarily 
exists within this frequency band. Backward waves in a lat- 
tice correspond to negative dispersion, i.e., the group velocity 
(the derivative of the eigenfrequency with respect to the wave 
number) is in the opposite direction of the phase velocity. Neg- 
ative dispersion for a lattice is quite common in high frequency 
bands (kd > it, where k = 2tt/X is the wave number in the 
host medium). However, at low frequencies (kd < it) nega- 
tive dispersion is an abnormal phenomenon. When kd < it the 
wavelength A is high enough as compared to the lattice period 
and the SRR size and thus the homogenization of the lattice 
may be possible. In [^],[^] and [|l0|] the lattice at low frequen- 
cies was treated as a continuous medium and the concept of the 
negative dispersion is equivalent to the concept of negative ma- 
terial parameters. However, is it possible to describe the struc- 
ture formed by lattices of wires and SRR:s (studied in [||] and 
[||]) in terms of e and /1 within the resonant band of SRR:s? 
This question remains open since the analytical model [||] with 
which these material parameters were introduced is incomplete 
and can be wrong. 

In [j|] the negative refraction was observed only for waves 
whose electric field is parallel to the wires and the wave vector 
is perpendicular to the wires. Thus, the permittivity e consid- 
ered in [0J and [|[] are the xx— component of the permittivity 
tensor (we assume the wires are along the x-axis) and the per- 
meability fi considered there are the transversal component of 
the permeability tensor /i t = fi yy = \i zz (there are two or- 
thogonal SRR:s in each unit cell of the structure studied in [j^]). 
If one considers only the propagation in the transversal plane 
(i.e., the y — z plane in our case), one can neglect the spatial 
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dispersion in the wire lattice and consider it as a medium with 
negative permittivity [g|. In the case of perfectly conducting 
wires one has [^] : 

e = l--f (1) 

where ui p is the analogue of the plasma frequency. In [Q, 
Q, [|l2j] the following frequency dependencies of the effec- 
tive medium parameters were suggested for their structure when 
both the SRR particles and the wires are made of real metal: 



e = 1 - 



,i = l- 



- JJUJ 



3iu 



(2) 



(3) 



where uj mp is the analogue of the resonant frequency of a mag- 
netic plasma, u e is the electronic resonant frequency and ui m is 
the magnetic resonant frequency (i.e.,, the resonant frequency 
of the magnetic polarizaibility for an individual SRR). These 
formulas correspond to the conventional quasi-static model for 
the permittivity and permeability, which takes into account only 
the local interaction in the structure under the homogenization 
(i.e., the electromagnetic interaction between the SSR:s and 
wires are neglected). In the present paper, we take into account 
the electromagnetic interaction (including the plane-wave in- 
teraction) between the SSR:s and the wires. We consider the 
lossless case (corresponding to 7 = 0, u> e = in Pendl and 
Pend2). 

In the present paper, we also discuss the following two ques- 
tions : 

• Is it possible to consider the structure formed by lattices 
of wires and SRR:s as a medium that its e is the same as 
the effective permittivity of a wire medium and its /j, is the 
same as the effective permeability of the lattice of SRR:s 
(i.e., is it possible to neglect the electromagnetic interac- 
tion between the SRR:s and the wires when calculating the 
material parameters of the whole structure)? 

• How to find the frequency region within the resonant band 
of SRR scatterers in which the homogenization of the 
whole structure is allowed? At some frequencies within 
the SRR resonant band the homogenization model may 
give so high values of the effective permeability that the 
product |e/i| becomes very large. Though the wavelength 
A in the host medium is much larger than the lattice period 
d, the effective wavelength A/ \f\ep\ of the eigenmode at 
these frequencies can be of the order of d or even smaller 
and thus the homogenization model becomes contradic- 
tory. However, within the band of the resonance there are 
frequencies at which the product |e/x| is not so large. Thus , 
the problem is how to separate quantitatively the frequen- 
cies at which the homogenization is forbidden from the 
frequencies at which the homogenization is allowed. 

The first question has been considered briefly in [ p"3| ] . It was 
indicated that the structure considered in [Q] was fortunately 
built so that each SRR is located exactly at the center of two 
adjacent wires and thus there is no quasi-static interaction be- 
tween the wire lattice and the SRR lattice (i.e., the magnetic 



fields produced by the two adjacent equivalent line currents can- 
cel out at the center where the SRR is located). If one considers 
e and ^ as quasi-static parameters (as was done in [[Ojl), the ab- 
sence of the quasi-static interaction should lead to the following 
result: the effective permittivity of the structure is identical to 
the effective permittivity of the lattice of wires and the effective 
permeability of the structure is identical to the effective perme- 
ability of the lattice of SRR:s. However, we will show in the 
present paper that this is not correct since the electromagnetic 
interaction between the wires and the SRR:s in such a structure 
is not quasi-static (or local) and will dramatically influence the 
effective permittivity. The consideration of structure as lattice 
of wires positioned in a host negative magnetic suggested in 
[ |l6| l is not adequate, because it also does not describe electro- 
magnetic interaction correctly. 

For the second question, one must be very careful in the ho- 
mogenization of the complex structure studied in [0], [^|] and 
[|To||. In fact, the results of [||] can not be interpreted quantita- 
tively in terms of the permittivity and permeability used in [Q] 
and [g], [g], [[§ and [0. This has been revealed in [Q. 

In the present work we develop an analytical model for a 
structure similar to the one studied in ^ (i.e., formed by com- 
bined lattices of infinitely long wires and split-ring resonators). 
The model allows the structure to be homogenized and its valid 
frequency domain to be identified. A self-consistent disper- 
sion equation is derived and then used to calculate correctly the 
effective permittivity and permeability in the frequency band 
where the lattice can be homogenized. The backward wave re- 
gion (which is a part of the resonant band of SRR) in which 
both the effective permittivity and permeability are negative is 
analysized. Since our homogenization is based on the disper- 
sion characteristics of the structure, the extracted material pa- 
rameters are not quasi-static (i.e., not only the local but also 
the non-local interactions are considered among the elements 
of the structure). Our results have shown that the homogeniza- 
tion is allowed only over part of the resonant band of the SRR 
scatterers and the homogenization is forbidden in a sub-band 
inside the SRR resonant band (even though the frequencies of 
this sub-band are low and the spatial dispersion exists there). 

II. SRR WITH IDENTICAL RINGS 

The SRR particle considered in [^J and [Q] is a pair of two 
coplanar broken rings. Since the two loops are not identical the 
analytical model for this particle is rather cumbersome (SRR 
models more complete than those suggested in Q| and [Q| have 
been developed in [ |l5| ] and Jl^]). It is not correct that the SRR 
particle can be described simply as a resonant magnetic scat- 
terer [ |l9| ] . The structure considered in [j|] (if homogenized) has 
to be described through three material parameters: e, fj, and n 
(the magnetoelectric coupling parameter). It was shown in [ |l9| l 
that a SRR is actually a bianisotropic particle and the role of 
bianisotropy is destructive for negative refraction. 

The bianisotropy is not the only disadvantage of this SRR 
particle. Another disadvantage is its resonant electric polar- 
izability (i.e., the polarization produced by the electric field), 
which was also ignored in and [H], Electric polarizability 
resonance occurs at frequencies very close to the resonant fre- 
quency of the magnetic polarizability [|l5|]. If the SRR:s are 
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Fig. 1 

Left: coplanar SRR. Right: parallel SRR. The concept of the 

MUTUAL CAPACITANCE IS ILLUSTRATED BY INDICATING THE CHARGES IN 
THE POLARIZED RINGS. 



made of real metal, the resonant electric polarizability will lead 
to a dramatic increase in the resistive loss. The resonant electric 
polarizability also makes the analytical modelling of the whole 
structure very complicated. 

A modified SRR which does not possess bianisotropy was 
proposed in [Jl9|], This SRR also consists of two loops but they 
are identical and parallel to each other (located on both sides of 
a dielectric plate in practice). Fig. |l| shows two kinds of SRR. 
The left one is the SRR considered in and [Q], The right 
is the SRR introduced in [ |l9l | and the one considered in the 
present paper. It has been mentioned in [|l9|] that the magnetic 
resonant frequency of their SRR is lower than that of the SRR 
considered in [|^] and [|J (for the same size). This is because the 
mutual capacitance C mu t between the two parallel broken rings 
is now the capacitance of a conventional parallel-plate capacitor 
and is significantly higher than the mutual capacitance of two 
coplanar split rings considered in [Q]. This fact is illustrated 
in Fig. [j]: if the upper half of ring 1 is charged positively the 
negative induced charges appear in the upper half of ring 2. 
The same situation happens for the SRR shown in the left part 
of this figure, but not so effectively since the strips are coplanar 
and weakly interacted. Therefore, the homogenization is more 
appropriate within the resonant band of the particle since the 
ratio of the particle size to A is small. 

There is one more advantage that was not mentioned in [ |l9| ] 
for the SRR of two parallel broken rings. The resonances of 
the electric and magnetic polarizabilities of this particle do not 
overlap in frequency. Actually, particle suggested in [19| is a 
special case of the bi-helix particle introduced (with the aim to 
create novel low -reflective shields) and studied in [[l7|]. Unlike 
the SRR considered in [|l9|| this bi-helix particle contains four 
stems orthogonal to the loop planes. Note that the theory of [ |l7| ] 
remains valid even if the length of these stems becomes zero. 

Both rings 1 and 2 (see the right diagram of Fig. |lj) have the 
same radius r and area S = ■ 



The impedance Z§ for each 



broken ring can be calculated by 



Z 



1 



where Y^ng is the admittance for the corresponding closed ring 
and y sp iit is the admittance for the split (associated with the 
capacitance between the two broken ends). The magnetization 
arises due to the magnetic field orthogonal to the ring plane 
(i.e., the xy plane in Fig. [l]). Resonant electric polarization is 
caused by the y— component of the external local electric field 
E loc . The rc-component of E loc has no influence over the res- 
onant polarization and can be neglected jl8|]. Thus, voltages 
(electromotive forces) £ H and £ E will be induced in each loop 
by the external local electric and magnetic fields, respectively 
L7l], [tTF 



and 



, E _ 4rJ{(fcr)Z 



7"/(Xi'mg + ^split) 



(4) 



(5) 



where J{ is the derivative (with respect to the argument ) of 
the Bessel function, and A\ is one of the so-called King's coef- 
ficients known in the theory of loop antennas (see e.g. Jl8||). 
Then the induced currents 7x,2 (due to the changing of the 
charges at the tips of the split arms) at the split gaps of rings 
1 and 2 satisfy the following equations, [pTQ: 



IlZ mu t 



= £ E 
= £ E 



■£ H , 



(6) 
(7) 



Yri 



split 



where Z mut is the mutual impedance of the two broken rings. 

It is clear from ppp and mmm that there are two eigenmodes 
of currents in the SRR. The first mode corresponds to I\ = I2 
when the electric dipole moment of the SRR is zero and the 
magnetic dipole moment m = mzo is twice the magnetic mo- 
ment of a single ring (see Fig. [j]). The second mode corre- 
sponds to 1\ = — 12 when the magnetic moment is zero and the 
electric dipole moment p = pyo of the SRR is twice the elec- 
tric dipole moment of a single ring. The first mode in the SRR 
is excited by the local magnetic field, and the second mode is 
exceited by the local electric field. The more the two rings are 
mutually coupled, the more is the difference between the reso- 
nant frequencies of the electric and magnetic resonantors [ 17 1. 
The electric resonant frequency (at which the electric polariz- 
ability resonates) is always higher than the magnetic resonant 
frequency [[L7|]. The relative difference of these two resonant 
frequencies may exceed 50% [p"7|]. If the distance h between 
the two parallel broken rings is very small, we can assume that 
the mutual coupling is so strong that the electric polarizabil- 
ity of the SSR is negligible within the frequency band of the 
magnetic resonance. Note that h is also the thickness of the 
dielectric layer between the two parallel broken rings in Fig. 
[j] since the rings are assumed to be perfectly conducting with 
infinitesimal thickness. 

Therefore, unlike the SSR considered in [f| and [[§, the SSR 
suggested in [ |l9| ] is appropriate for creating artificial magnetic 
resonance without resonant electric properties within the fre- 
quency band of interest. This is the reason why we choose the 
SSR shown in the right part of Fig. [l] to study in the present 
paper (the analysis will be much simpler). 

In the present paper we introduce an analytical model for an 
individual SRR particle used in JT^]. This model is simpler than 



4 



IEEE TRANSACTIONS ON ANTENNAS AND PROPAGATION, SPECILA ISSUE ON METAMATERIALS, 2003 



the one considered in [ |17| ] (due to the absence of the stems). 
Assume that the SRR shown in Figure [l] is made of perfectly 
conducting strip and is excited by magnetic field H l ° c . Also 
assume that the dielectric plate separating the two parallel rings 
has the same permittivity as the background medium (then we 
can avoid the influence of the dielectric plate which can be very 
strong if there is a mismatch in the permittivity). The model 
used in this section is quasi-static since it refers to an isolated 
particle of small size (with respect to the wavelength). 

If the non-uniformity of the azimuthal current distribution in 
both rings and SRR:s can be neglected, the magnetic polariz- 
ability can be written as 



2I> S 



H 



loc 



H 



loc 



where I = I\ = 1%. From ppp and mag it follows that 

£ H -ju>Sn H l ° c 



Zq + z n 



Zq + 2 mu t 



(8) 



(9) 



Calculating the total impedance of the loop by taking into 
account the mutual coupling of the loops (as it was done in [ |l5| ] 
and Jl^ ] for SRR of coplanar rings), we obtain 



Ztot — Zo + Z mut — ju>(L + L mut ) + 



R r 



Here R r is the radiation resistance of the whole particle, L is 
the ring inductance: 

32r 



L = //or log 



where w is the width of the strip from which the ring is made of, 
and L mut is the mutual inductance of the two parallel coaxial 
rings: 

3£ 2 \, 4 



1 



4 y log c 



where £ = h/2r. The total capacitance Ctot attributed to the 
split can be calculated (taking into account the capacitive mu- 
tual coupling; cf. [ |T5| ] and jl9tl) as half of the mutual capaci- 
tance formed by the two parallel rings 



Cu 



C mut _ e ewnr 
2 ~ 2h 



In this formula the capacitance of the split is neglected since it 
is small as compared to C mu t. 
From def and tok we obtain 



2/igS 2 



(L + L mut )(§-l)-j% 



where 



(10) 



(11) 



Ctot(L + Lmut) 

In a similar way one can show that the electric polarizability 
resonates at the frequency lj\ (see also [|l7|]): 

Ctot{L — imut) 



b 



a : 



+- +■ 



y Wires SRR 



g =1^° H^o H»e H^ 
4 4 + 4 + 
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Fig. 2 

Left: front view of the lattice of SRR:s (shown as disks; their 
magnetic dipoles are indicated with arrows) and straight 
wires . Right: top view. The wave propagation is in the y — z 

PLANE. 



and o^o < uj\. 

We will also use the following result of amm: 



Re 



1 



(L + £ rout )(£i-i) 



2^ 2 



(12) 



The radiation resistance R r can be found from the following 
condition [El, & [pi: 



Im 



1 



(13) 



In the dispersion equation for a lattice, R r cancels out and does 
not influence the result. 

III. The structure 

The structure we study in the present paper is similar to 
the one studied experimentally in Hi, however, instead of the 
coplanar SRR:s we use the parallel SRR:s (as described in the 
previous section). 

When the wave propagates along the z axis the electric field 



(n y , n z ) (for the reference wire we have 



0). The 



magnetic field excites those SRR:s which are parallel to the 
x — z plane. Their magnetic moments are parallel to the y 
axis. Then the lattice can be considered as a set of 2D grids 
parallel to the x — y plane and orthogonal to the propagation di- 
rection. Each grid contains magnetic and electric polarizations. 
The magnetic moments as well as the currents are tangential 
to the grid plane, and each grid can be considered as a sheet 
of surface magnetic moment M = Myo and a surface electric 
current J = Jxo (or surface electric polarization P = J / juS). 
Similar situation holds for the wave propagation along the y 
axis. Then one has M = Mzq and J = Jxo, and the elec- 
tric and magnetic polarizations for each 2D grid (parallel to the 
x — z plane) is again tangential to the grid and orthogonal to 
the propagation direction. The wire lattice and the SRR lattice 
have the same periods along the y and z axes and are denoted 
as b and d, respectively. The period of the SRR lattice along the 
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x axis is denoted as a. Fig. [j] shows the two orthogonal sets of 
SRR:s separated with each other by a/2 along the x axis. This 
separation plays no role in our model since we do not consider 
the electromagnetic interaction between these two orthogonal 
sets of magnetic dipoles. When the wave propagates along the 
z axis or y axis, one of the two sets of SRR:s is not excited and 
the interaction is completely absent. 

In the case b = d such a structure behaves (within the fre- 
quency band where the homogenization is possible) like a uni- 
axial magneto-dielectric medium with relative axial permittiv- 
ity e xx (mainly due to the presence of the wires) and relative 
transversal permeability fi yy = fj, zz = /j t (mainly due to the 
SRR particles). This indicates that in order to find the effective 
material parameters of the whole structure we can consider only 
the case of the normal propagation (along the z or y axis). 

Note that the transversal permittivity and the axial perme- 
ability of the structure are equal to those of the background 
medium. For simplicity we assume this background medium 
is vacuum. 

We will see that the SRR:s strongly interact with the wires 
at each frequency. Their interaction is not quasi-static and in- 
fluences the propagation constant starting from zero frequency. 
In this way it influences the material parameters of the whole 
structure. 

IV. Dispersion equation 

Let the wave propagate along the z axis with propagation 
factor (3 (to be determined). Consider the whole structure as a 
set of parallel 2D grids which are parallel to the x — y plane and 
denote the surface magnetic moments M(n z ) and the surface 
currents J(n z ) at those grids numbered n z . Then we choose an 
arbitrary SRR in the grid with n z = as the reference particle 
and an arbitrarily chosen wire (in the same grid) as the reference 
wire. 

When we evaluate the magnetic moment m of the reference 
SRR (which is related to the surface magnetic moment M by 
rn = Mab), we take into account its electromagnetic interac- 
tion with all the other SRR:s following the work of [ p4| ] where 
a simple model of 3D dipole lattice was suggested. As to the 
influence of the wires to the reference SRR, we can replace 
each grid of wires with a sheet of current J(n z ) because of 
the absence of the a quasi-static interaction between SRR:s and 
wires. This gives the well-known plane-wave approximation 
of the electromagnetic interaction in lattices (see [ p4| ] and the 
references cited there). 

When we evaluate the current 7 of the reference wire (which 
is related to the surface electric current by 7 = Jb), we take 
into account its interaction with all the other wires following 



the work of [ 20 ] where a simple model of doubly-periodic wire 
lattice was suggested. The influence of the SRR lattice on the 
reference wire is taken into account under the plane-wave ap- 
proximation and the reciprocity principle is satisfied. 

Each sheet of electric or magnetic polarization produces a 
plane wave [Bin. Since M(n z ) and Jin z ) satisfy 



M(n z ) 
J(n z ) 



(14) 
(15) 



we can write the following relations for the x— component 
of the electric field (produced by all the sheets of mag- 
netic moment M and acting on the reference wire) and the 
y— component of the magnetic field (produced by all the sheets 
of current J and acting on the reference SRR): 



E? = 



H, = 



E 

n z — — oo 
oo 



J 



(16) 



^ sign(n z )^e- jn ^ d - jk \ n ^ d (17) 

n z — — oo ~ 



Both series can be analytically carried out and we easily obtain 

uj M sin (3d 



= -- 



77: 



2 cos kd — cos /3d 
jj sin/3d 
2 cos kd — cos (3d 



(18) 
(19) 



The local electric field acting on the reference wire is the sum 
of El 1 and the contribution of the wires: 



E 



loc 



E 



M 



C W I = E 



A I 



C w bJ 



(20) 



where C w (the interaction factor of the wire lattice) was deter- 
mined in [EG] : 



2b 



sin kd kb ( kb \ kb 
— — + — log — + 7 + j — 

cos kd — cos (3d it \ in J 2 

(21) 

Here 77 is the wave impedance of the host material and 7 = 
0.5772 is the Euler constant. 

The local magnetic field acting on the reference SRR is the 
sum of Hy 1 and the contribution of the SRR particles: 



rrloc ttJ 

H y -Hy- 



dra = Hi 



C d abM 



(22) 



where d (the interaction factor of the lattice of magnetic 
dipoles) is given by [E4fl 



C d = 



ujq 



+ 3 



i kd 



2abr] &nmuQ 2abr\ cos kd — cos (3d 



(23) 



A similar relation has been given in [ p4[ | for a lattice of electric 
dipoles (the only difference as compared to bbb is the factor r/ 2 ). 
Here q denotes the real part of the dimensionless interaction 
factor of a 2D grid of dipoles with periods a, b. In [ p4| ] the 
closed-form expression for q is given for the case a — b: 



1 

90 = 2 



kas 



— sin kc 



where the number s is approximately equal to 1/1.4380 = 
0.6954. Relation bbb is very accurate for the case d ^> a, and 
in the case d — a its error is still quite small pi}]. 

The responses of the reference SRR and the reference wire 
to the local fields can be written as 



ra 
I 



,H 



loc 



E 



loc 



(24) 
(25) 
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where Z w is given by 

_ krj 



(26) 



where ro is the effective radius of the wire (7*0 = w/4 if made 
from a strip with width w). 

To obtain the dispersion equation we substitute formulas 
eloc, hloc, aaa, bbb, esum and hsum into fefm and fefw. Since 
m = Mab and I = Jbwe obtain the following system of equa- 
tions: 



M 



1 



A 



3: 



a mm 67T/xo 



.BJ 
^ ab 
= -ujBM 



(27) 
(28) 



where 



and 



A = 



sin kd 



2abr\ \ cos fed — cos /3d 
1 sin /3d 

i3 



2 cos fed — cos /3d 

The parameter £> describes the interaction between the currents 
in the wires and the magnetic moments of the SRR:s. B is not 
a quasi-static parameter even at low frequencies since it does 
not approach zero at zero frequency. Its presence in the disper- 
sion equation strongly influences the result for the propagation 
constant /3 at all frequencies. 

Relations ima and zwire lead to the cancellation of the imag- 
inary part at the left-hand side of fel and the real part at the 
left-hand side of fe2. Thus, system of equations fel and fe2 
gives the following real-valued dispersion equation: 



(Im(Z w ) - \m{C w )) 



Re 



1 



-A 



lu 2 B 2 
ab 



It can be re-written as the following quadratic equation with 
respect to cos (3d: 



2a6. 

LOT] 



-Re 



( 1 ) 


- q 







(cos kd — cos /3d) sin kd 



kb, b 
1 + — log 

ir 27rro 



1 2 fed - 



cos 2 /3d + 1 = (29) 



There are two roots j3i$ for the dispersion equation exp at each 
frequency. One of them is exactly equal to fc (the wavenumber 
in the host medium). This root corresponds to the wave with 
polarization E = E y yo and H = H x xq. This wave excites 
neither wires nor SRR:s and does not interact with the structure. 
Another root corresponds to the wave with polarization E = 
E x xq and H = H y yo- This is the interacting wave which is 
of interest. In our dispersion curves we keep both solutions of 
exp. 

V. Dispersion curves 

As numerical examples we choose the following parameters 
for the structure shown in Fig. ||: the size of SRR particle 
(outer diameter of the rings) is D = 3.8 mm, the width of 
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Fig. 3 

DISPERSION PLOT OF THE STRUCTURE FOR d = 2a = 2b = 16 MM (THICK 

LINES). THIN lines coincide with the dispersion plot of host 

MEDIUM AND CORRESPOND TO THE WAVE POLARIZATION 

E = E y y ,H = H x x . 



the strip (forming the rings) is w = 1 mm, the radius of wire 
cross section is ro = 0.2 mm, the distance between the rings 
(which is chosen so that the resonance of a mm is at 6 GHz) is 
h = 0.84 mm. Lattice periods a = b = 8 mm and d = 16 
mm are chosen in our first example. In the second example, we 
choose a = b = d = 8 mm. In the third example, we choose 
a = b = d = 4 mm. 

Fig. U gives the dispersion curve in a form commonly used in 
the literature of photonic crystals (see e.g. [p6[]). It represents 
the dependence of the eigenfrequencies on the normalized prop- 
agation factor fid/ tt over the first Brillouin zone (0 < {3 < ir/d) 
for the case when d = 2a = 2b = 16 mm. Straight lines corre- 
spond to non-interacting waves. Curved lines correspond to in- 
teracting waves. The only difference of this curve as compared 
to the well-known plot for the wire medium (see e.g. and 
1^5]]) is the mini -band at about 6 GHz, in which the group veloc- 
ity is negative. In this narrow frequency band wave propagation 
is prohibited in the lattice of wires. Therefore, the mini-band is 
due to the presence of SRR:s and the resonant magnetization of 
the SRR lattice. The eigenfrequencies (associated with waves 
with E = E x xq and H = H y yo) within this pass-band are 
shown by the crosses. 

The resonant pass-band becomes wider if the period d de- 
creases. From Fig. [f]one can see the dependence of the propa- 
gation factor on the frequency in the vicinity of the SRR reso- 
nance for the case a = b = d = 8 mm. The solid line in this 
figure corresponds to the non-interacting wave. 

The backward-wave region corresponds to the frequencies 
5.980-6.045 GHz, whereas the resonant frequency of Re(a mm ) 
given by omega is 6.000 GHz. Re(a mm ) becomes negative at 
6.000 GHz. Thus, within the backward-wave band Re(a mm ) is 
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Normalized propagation factor of the backward wave versus 

FREQUENCY FOR d = a = b = 8 MM. 



mainly negative. The group velocity of the backward wave is 
relatively small (it approximately equals 5.3 • 10~ 3 c, where c is 
the speed of light). 

In order to understand whether it is possible to homogenize 
the structure at the frequencies when the backward-wave region 
exists, we studied (3 for both propagating and decaying modes. 

Fig. H shows the frequency dependence of both the real and 
imaginary parts of the normalized propagation factor (3d/ V for 
the case a = b = d = 8 mm. It is clear that outside the SRR 
resonant band the eigenmodes of the structure are the same as 
those of the wire medium The thick straight line corre- 
sponds to the non-interacting wave. 

From Fig. || one can see that the eigenmodes within the fre- 
quency band 5.92-5.98 GHz are complex. The lower limit of 
the backward wave region (5.980-6.045 GHz) is the upper limit 
of the complex-mode band. Complex modes cannot exist in 
continuous media. These modes are known for electromagnetic 
crystals with different geometries (see, e.g. [p0|], [p3[|). These 
are decaying modes though the real part of the propagation fac- 
tor is Re/3 = ir/d. The existence of this real part of (3 reflects 
the fact that the directions of the currents in the wires are alter- 
nating along the propagation axis (two adjacent currents have 
opposite directions and this can be interpreted as the phase shift 
7r between them due to the real part of the complex propagation 
factor). 

Therefore, the homogenization is possible within one (the 
upper) half of the SRR resonant band but impossible within an- 
other (the lower) half of the SRR resonant band (though for 
these frequencies the structure periods are much smaller than 
the wavelength in the background medium). 

VI. Homogenization 

Let us try to consider the structure (in the case a = b = d) 
as a uniaxial magneto-dielectric medium. Then the interact- 
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normalized propagation factor versus frequency for 
d = a = b = 8 mm. Thick line corresponds to the propagation of 
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ing wave (propagating along the z axis with E = E x Xq,H 
H y yo) also satisfies the following constitutive equations: 



I), = <,,/:, + /'; 



bulk 



By — ^qH v + M y — noiitHy 



Define the following ratio: 

^>bulk 



a = 1] 



1 {ex 



I) Er. 



Af bulk 77 (fi t - 1) Hy 



(30) 



where E x and H y are the field components averaged over the 

cubic cell a x a x a. Pi? ulk and M? ulk are the bulk electric and 

* y 

magnetic polarizations related with the surface current J and 
surface magnetic polarization: 



^jbulk 



j 

jujd 



/.xbulk _ M 

y ~ d 



From Maxwell's equations we easily obtain 

E x 

Hy 

Substituting aaa and zwire into fe2, we obtain: 
■qj rjB 




(31) 



a(u>) 



jujM j(Z w - C w ) 
it sin (3d 



7r sin kd + kb log (cos kd — cos (3d) 



(32) 



27rr 



8 



IEEE TRANSACTIONS ON ANTENNAS AND PROPAGATION, SPECILA ISSUE ON METAMATERIALS, 2003 



From fat and imped it follows that 



1 5 r 



a 



Oh - 1) V e 



(33) 



In the above equation P(lo) is already known from the disper- 
sion curve and k = uj^/eofio = w/c. 

Equating the propagation factor j3 to the value to-y/eofioe xx fit> 
we obtain 

Substituting this expression for fit into alpha, we obtain: 



c (w) 



_c§_ 



ka(tu) 



(35) 



After e xx 
rhu) 



is found, we then evaluate fi t through (cf. eps and 



wq(lj) 



1 



(36) 



We have taken into account the non-local interaction in the 
structure in formulas eps and mumu though the effective per- 
mittivity and permeability are introduced as the parameters re- 
lating B. D with E, H at the same point. Therefore, unlike for- 
mulas Pendl and Pend2, our material parameters are not quasi- 
static. 

Frequency dependencies of both e xx and fi t are shown in Fig. 
|6| for the case a = b = d = 8 mm. Outside the resonant band of 
the SRR particles the frequency dependence of the permittivity 
repeats the known result for wire media (treated as artificial di- 
electric media) The permeability is practically equal to the 
unity outside the resonant band of SRR:s. Within the complex- 
mode band the homogenization is forbidden [53], J23| ] and this 
frequency region is removed in this figure (both e xx and fit cal- 
culated through eps and mu are complex within this band). 

Let us consider the resonant frequency behavior of the ma- 
terial parameters in details. Fig. [7] shows the same curves as 
those in Fig. however, in another scale starting from the up- 
per limit of the complex-mode band. From Fig. [7] one sees 
that the permittivity and permeability are both negative between 
5.970 and 6.020 GHz. From Fig. §it follows that the backward 
wave propagates between 5.980 GHz and 6.045 GHz. Thus, the 
backward-wave region almost coincides with the region where 
both the permittivity and permeability are negative. Note that 
this coincidence is only approximate. 

Also we have indicated in Fig. [7| the point at which per- 
mittivity and permeability are equal (at about 6.005 GHz). At 
this frequency, the medium is impedance-matched with the free 
space (this is useful for some applications) and the values of e xx 
and fi t are not very high (the homogenization is then allowed). 

As a main result, one can see from Fig. [7] that the permit- 
tivity does not follow (within the resonant band) the law (even 
qualitatively) suggested in [||], [^] and [|I(J|. it's the frequency 
dependence of the permittivity is non-monotonous (from Fig. ^ 
one sees that e xx decreases over 5.96-6.09 and increases after 
6.090 GHz as the frequency increases). 
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PERMEABILITY (SOLID LINE) FOR THE CASE a = b = d = 8 MM. 



In the theory of continuous media one can prove that both 
the permittivity and permeability must grow as the frequency 
increases in the lossless case ||2^]. In our case the permeabil- 
ity grows everywhere as the frequency increases (until the first 
spatial resonance of the lattice, i.e., kd = it when it loses the 
physical meaning). Thus, the frequency behavior of the per- 
meability is normal. However, the permittivity grows as the 
frequency increases only at the frequencies when the magne- 
tization of SRR is small and the interaction of the SRR:s and 
wires is negligible. Within the band of the backward wave the 
permittivity decreases as the frequency increases. Therefore, 
the homogenization procedure we have developed is not com- 
pletely consistent with the theory of J27[]. The reason for this 
disagreement is that the material parameters considered in [ ^7| [ 
are quasi-static (i.e., the polarization of the medium at a given 
point is determined by the field at this point) while our model 
takes into account the non-local interaction of the SRR lattice 
and the wire lattice. We found that the visible difference be- 
tween the quasi-static model and our model is within the SRR 
resonant band. However, the influence of the non-local interac- 
tion is revealed in the permittivity of the wire lattice disturbed 
by the presence of the SRR:s. 

The lattice of infinite wires is spatially dispersive at all fre- 
quencies since the wires are longer than any possible wave- 
length. When the wave propagates strictly in the plane or- 
thogonal to the axis of the wires one can still neglect the spa- 
tial dispersion since all parameters are independent of the x- 
coordinate. Thus, the problem is two-dimensional and possible 
to be homogenized |^]. However, if there is a lattice of scatter- 
ers with which the wires interact, the situation becomes quite 
different (even for propagation orthogonal to the wires). Here 
the problem is not two-dimensional and the wire current is in- 
fluenced by all the SRR particles positioned along its infinite 
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length. It results in the abnormal frequency behavior of the ef- 
fective permittivity of the structure. 

VII. Conclusion 

In the present paper we have developed an analytical model 
for a structure similar to the one for which the negative re- 
fraction at microwave frequencies was first observed (formed 
by combined lattices of infinitely long wires and split-ring res- 
onators) [fj]. We have derived a self-consistent dispersion equa- 
tion and studied the dispersion properties of the lattice. The 
explicit dispersion equation clearly confirms the existence of 
the narrow pass-band within the resonant band of the split-ring 
resonators. In this pass-band the group and phase velocities of 
the propagating wave are in opposite directions (i.e., a back- 
ward wave). Negative refraction is the direct consequence of 
the backward wave (instead of the negative values of permittiv- 
ity and permeability as thought in many works [|l4|]). Thus, the 
experimental observations of are consistent to our model. 

Our model also allows the homogenization of the compos- 
ite structure. The obtained dispersion curves have been used to 
calculate correctly the effective permittivity and permeability in 
the frequency band where the structure can be homogenized. It 
is interesting to see that the dispersion curves agree well with 
the Veselago theory which predicts backward waves when both 
permittivity and permeability are negative. Outside the reso- 
nant band of the SRR particles, the effective permittivity of the 
whole structure is the same as that of the wire lattice and the 
effective permeability is equal to 1 . However, within the SRR 
resonant band, there is a sub-band where the homogenization is 
forbidden since the complex mode satisfies the dispersion equa- 
tion at these frequencies. We found that the frequency region in 
which both e and /i are negative coincides approximately with 



the backward wave band. In this region the frequency depen- 
dence of the effective permittivity is abnormal. We interpret 
this as the result of the low-frequency spatial dispersion which 
is inherent for the wire medium in the presence of the resonant 
scatterers. 
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